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Abstract 

The diverse functionality emerging at oxide interfaces calls for a fundamental understanding of 
the mechanisms and control parameters of electronic reconstructions. Here, we explore the evolu¬ 
tion of electronic phases in (LaA103)M/(SrTi03)Ar(001) superlattices as a function of strain and 
confinement of the SrTi03 quantum well. Density functional theory calculations including a Hub¬ 
bard U term reveal a charge ordered Ti^+ and Tp+ state for = 2 with an unanticipated orbital 
reconstruction, displaying alternating dxz and dyz character at the Ti^+ sites, unlike the previously 
reported dxy state, obtained only for reduced c-parameter at asxo- At aLAO c-compression leads 
to a Dimer-Mott insulator with alternating dxz^ dyz sites and an almost zero band gap. Beyond 
a critical thickness of = 3 (asTo) and A^ = 4 (aLAo) an insulator-to-metal transition takes 
place, where the extra e/2 electron at the interface is redistributed throughout the STO slab with 
a dxy interface orbital occupation and a mixed dxz + dyz occupation in the inner layers. Chemical 
variation of the SrTi03 counterpart (LaA103 vs. NdGa03) proves that the significant octahedral 
tilts and distortions in the SrTi03 quantum well are induced primarily by the electrostatic doping 
at the polar interface and not by variation of the SrTi03 counterpart. 
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Oxide interfaces show surprisingly rich electronic behavior that is not present in the bulk 
compounds. Among the different systems, the (OOl)-oriented LaAlOs/SrTiOs interface has 
been most prominent, displaying two-dimensional conductivity, superconductivity and mag¬ 
netism, as well as confinement-induced metal-to-insulator transitions^. The majority of 
experimental studies have explored thin LaAlOa (LAO)-films on SrTiO 3 ( 001 ). In contrast, 
superlattices containing two n-type interfaces with a LaO and Ti02 stacking have been in 
the focus of theoretical studieS^^. The latter predict that the excess charge of 0.5e at the 
polar interface results either in a charge ordered phase with d^y orbital polarization or a quasi 
2DEG where the excess charge of 0.5e is delocalized throughout the SrTiOs (STO) quantum 
well. While the experimental realization of this type of superlattice remains a challenge as 
it requires precise control of termination at both interfaces, there is renewed interest in such 
systems in view of confinement effects in electrostatically doped oxide quantum wells. To 
this end, recent experimental studies on related GdTiOs/SrTiOs super lattices by Moetakef 
et al^ reported the formation a high density 2DEG where the carrier concentration can 
be modulated by the thickness of the STO quantum well. Moreover, a metal-to-insulator 
transition and ferromagnetism is induced for STO thickness below two unit cells (containing 
two SrO layers)!^. In contrast, the analogous SmTiOa/SrTiOa superlattice (SL) remains 
metallic down to a single SrO layei®. This distinct behavior of the two systems has been 
related to differences in octahedral tilts and distortions, induced by GdTiOs or SmTiOs. 
These findings have motivated us to revisit the LAO/STO superlattices and explore sys¬ 
tematically the types of electronic and orbital reconstructions that can be realized as a 
function of STO confinement. We note that LAO/STO SLs have the advantage that the 
transition metal ion is only confined to the STO part, therefore emerging effects, includ¬ 
ing magnetism can be unambiguously associated with the STO quantum well (QW). We 
relate the electronic behavior to the structural changes by fully considering octahedral tilts 
and distortions as well as the role of the vertical lattice parameter. We note that pre¬ 
vious work® has addressed predominantly tetragonally distorted superlattices, where due 
to symmetry restrictions, only a structural optimization of the position of the ions was 
performed. GdEeOs-type distortions have been considered by Zhong and Kellj^^ only for a 
(LAO)m/(STO)jv( 001) superlattices with N = 3. Here, we have investigated the behavior 
as a function of A. In particular, a critical thickness for a metal-to-insulator transition in 
this system is determined. 
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Strain is a powerful tool to tune the electronic behavior and orbital polarization^HUSl^ 
so far only a few studies have addressed this effect in the LAO/STO system. Bark et al^ 
studied its influence on the properties of LaAlOs Aims on SrTiO3(001) by growing the system 
on substrates with different lattice parameters. Transport measurements show that under 
compressive strain the 2DEG is preserved, but with reduced carrier concentration compared 
to the system grown on SrTiO3(001), whereas for larger tensile strain the 2DEG is quenched. 
Furthermore, recent experiments have shown variations of carrier density and mobility of 
the 2DEG when exchanging LaA103 by other polar RBOs oxides (i?=La, Pr, or Nd, B=Al 
or Ga)P^. 
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FIG. 1: Side view of the relaxed structure, electron density distribution, integrated over occupied 
Ti 3d bands between Ep-l.S eV and Ep and band structure in (LAO)M/(STO)Ar (001) SL at asTO- 
Majority and minority bands are plotted in blue and orange, respectively. Character plotting 
highlights the occupied d band at the Ti^+ sites, (a) = 2 is a charge-ordered FM insulator with 

dxy orbital occupation at the Ti^+ sites for the c-parameter fixed to conserve the volume and (b) 
with dxz and dyz orbital occupation at Ti^+ at each interface layer for the optimized c-parameter. 
AE is the energy difference (in eV per Ti site) between the cases shown in (b) and (a);(c)-(d) 
For N — 3^4 STO layers, the charge-ordered state is suppressed and the system switches from 
insulating to conducting behavior with a dxy orbital at the Ti interface sites and dxz + dyz orbital 
occupation in the central STO layers. 

The goal of this theoretical study is to explore systematically the influence of strain, conflne- 
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ment and chemical variation on the electronic and orbital reconstruction in (OOl)-oriented 
(LAO)M/(STO)Ar snperlattices, where N is varied between 2-4. Octahedral tilts and distor¬ 
tions, considered to be an important factor in tuning the properties of the 2DEG at these 
oxide interfaces, were fully taken into account. Additionally, we also performed optimiza¬ 
tion of the ont-of-plane lattice parameter c. The effect of chemical control on the interface 
properties is addressed by replacing the polar LaAlOa with the polar NdGa03(NG0). 

(a) 2/4 (b) 2/4 (c) 3/3 (d) 4/2 



FIG. 2: Side view of the relaxed structure, electron density distribution, integrated over occupied 
Ti 3d bands between Ep-l.S eV and Ep, and band structure of (LAO)m/(STO )n (001) SL at apAO? 
N = 2 — 4. Majority and minority bands are plotted in blue and orange, respectively. Character 
plotting highlights the occupied d band at the Ti^+ sites, (a) For N = 2 a. Dimer Mott insulator is 
obtained for the strongly compressed c parameter with alternating dxz^ dyz orbital occupation at 
the Ti^-^+ sites; (b) for the optimized c a charge ordered state with dxz and dyz orbital occupation 
at the Ti^+ sites at each interface layer emerges; AE is the energy difference (in eV per Ti site) 
between the cases shown in (b) and (a); {c) N — 3 is similar to (b) with Ti^+ in the central layer; 
(d) an insulator to metal transition takes place for N — 4: with a dxy orbital polarization at the 
Ti interface sites and dxz + dyz orbital occupation in the inner STO layers. 
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Results 


Electronic and orbital reconstruction (LAO)m/(STO)7v( 001) 
SLs: influence of strain 

Results for asTO, corresponding to an underlying STO(OOl) substrate. 

Figure [^presents results for n-type (OOl)-oriented (LAO)M/(STO)Ar superlattices with vary¬ 
ing thickness of the STO quantum well, N = 2 — 4. Using the c parameter that corresponds 
to volume conservation of the corresponding bulk constituents, the SrTiOs bilayer {N = 2 
case) under tensle strain (osto) [Fig- [^a)] is a charge-ordered ferromagnetic insulator with 
alternating TU+ (0.59and TU+ (0.11/r^) sites with a d^y orbital occupation at the 
Ti^+ sites, consistent with previous workP®. Surprisingly, for the optimized c parameter of 
3.87 A (0.23 A larger) the charge-ordered pattern is preserved, however the orbital order at 
the Ti^"*" sites changes to d^z for TU'*' at the top interface and dy^ for TU"*" at the bottom 
interface [Fig. I3b)l- This is associated with a much narrower occupied Ti 3d band and 
an enhanced effective mass of 3.57 rUg, indicating an increase of electron correlation effects. 
Simultaneously the size of the band gap is enlarged from 0.16 eV to 0.42 eV. The magnetic 
moment of the TU'*' sites increases slightly from 0.59 fiB to 0.64/is. 

The mechanism of electronic reconstruction is distinct for the thicker STO quantum wells, 
N = 3 [Fig. [^c)] and N = 4 [Fig. Sd)]. Here, the extra e/2 charge from each interface 
is redistributed throughout the STO slab, resulting in an insulator-to-metal transition with 
a critical thickness Ng — 3. In both cases the charge-ordered state is suppressed. Simul¬ 
taneously, the orbital polarization is altered to predominantly d^y orbital character in the 
interface Ti-layer and mixed d^z + dyz in the central Ti layers, similar to previous findings 
in (LAO)m/(STO)jv superlattice tj^ * ^ * ^^ *^. 

Results for ai^AO, eorresponding to an underlying LAO(OOl) substrate. 

Figure shows the results for (OOl)-oriented (LAO)jvf/(STO)iv heterostructures with STO 
thickness ranging from A^ = 2 — 4. For the SrTiOs bilayer [N = 2) the ground state 
under compressive strain (olao) and optimized c parameter [see Fig. [^b)] is analogous to 
the one for usto [cf- Fig. Bb)]: a charge-ordered insulating phase with TF"'' and TF'*' in 
a checkerboard arrangement and a d^z, dyz orbital polarization at the Ti^"*" sites at each 
interface, albeit with slightly smaller magnetic moment (0.61 pb) at the Ti^"*" sites. This is 
accompanied by an increase of the band gap from 0.42 eV (osto) to 0.65 eV (olao)- 
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An interesting electronic state occurs for A = 2 upon compression of the c parameter 
of the superlattice to the bulk LAO value [see Fig. [^a)]: The charge-ordered state is 
quenched, leaving 0.5e at each Ti site (0.33/is) and alternating d^z, dyz orbital polarization 
on neighboring Ti sites in a checkerboard arrangement. Characteristic for this state is a 
decreased but nonzero band gap (0.07 eV). This so called Dimer Mott insulator (DMI), 
initially predicted within the framework of one-dimensional Hubbard ladderJ^^, was recently 
predicted for GTO/STO QWs for a single SrO layer in a GdTiOs matri:xP^. In fact the 
ground state for the optimized c-value can be regarded as emerging from the DMI phase 
by transferring 0.5e between each pair of Ti sites, leading to a TF"*" and TF"*" charge order 
with some residual d^z orbital occupation at the Ti^"*" sites and dyz at the Ti^"*" sites at one 
interface and the opposite orbital polarization in the other interface layer. 

To explore how the combination of exchange correlation functional (LDA vs GGA) and the 
Goulomb repulsion term influence the structural properties and electronic ground state, we 
have performed additional calculations for (LAO) 4 /(STO) 2 ( 001 ) at aLAO) using the LDA-I-C 
functional. Setting the lateral lattice parameter to the experimental bulk LAO value of 
3.79 A, the optimized c lattice parameter (3.83 A) is ~ 4% smaller than the GGA-l-7/ value 
(3.98A). The electronic ground state in Fig. [^b) with TF+ and TF+ in a checkerboard 
arrangement and a d^z, dyz orbital polarization at the TF+ sites emerges within LDA-l-7/ 
beyond U = 5.5 eV. Thus we can conclude that the electronic ground state is largely 
unaffected by the reduction of the c-lattice parameter within LDA-I-I7 and is stabilized for 
a slightly higher value of U than in the GGA-I-G calculation. 

Unlike the corresponding case at osto, the electronic and orbital reconstruction persists for 
N = 3 [see Fig. [^c)] i.e. a charge-ordered state with alternating d^z, dyz at the Ti^+ sites 
emerges at both interface layers, whereas the central layer becomes Ti^+. Compared to the 
case with N = 2 with optimized c, the band width is increased and the band gap is reduced 
to 0.14eV. The magnetic moment at the TU"*" sites is reduced to (0.54/is), indicating an 
enhanced delocalization. 

An insulator-to-metal transition takes place at a STO thickness N = 4 [see Fig. [^d)]. The 
suppression of the charge ordered state leads to a closing of the band gap. The extra e/2 
charge is delocalized throughout the whole STO quantum well with a preferential d^y orbital 
occupation at the interface layer and a d^z + dyz occupation in the central layers. 
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structural relaxations 


To investigate the relationship between electronic properties and strnctural distortions we 
have performed a quantitative analysis of the fully relaxed (LAO)jvf/(STO)iv ( 001 ) super¬ 
lattices. The B-O-B bond angles in the perovskite structure ABO 3 are closely related to 
octahedral tilts and determine orbital overlaps, and are therefore considered to be the con¬ 
trol parameter for metal-to-insulator transitions in bulk transition metal oxided^IES 
derstanding the origin of modifications of octahedral tilts and connectivity at interfaces in 
oxide superlattices is an emerging topic of considerable intereslP^. 



Oz-coordinate O z-coordinate 



Oz-coordinate Oz-coordinate 

FIG. 3: Quantitative analysis of B-O-B bond angle in (LAO)m/(STO)jv ( 001 ) superlattices for 
different STO thicknesses N (shaded region), (a) The bilayer N = 2 under tensile strain (osto)- 
(b) STO bilayer for compressive strain (aLAo)- (c)-(d) A = 3,4 STO layers. 


In Fig. 1^ we present the B-O-B bond angles of the superlattices as a function of the in¬ 
plane lattice constant asTO and ulao- The corresponding results for aNco are displayed 
in Fig. In the (LAO)m/(STO)jv ( 001 ) superlattices we identify strong reduction of 
the B-O-B angles from the bulk values. We remind that bulk SrTiOa has the ideal cubic 
perovskite structure {Pm3m) with B-O-B bond angles of 180°. Below 105 K it undergoes an 
antiferrodistortive phase transition to a tetragonal {lA/mcm) structur^, whereas LaAlOs 
attains a rhombohedral structure (i? 3 c) at all temperatures with B-O-B bond angles of 
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171.4°. The deviations of the B-O-B bond angles from the ideal 180° are most prominent 
in the STO quantum well. Moreover, there are pronounced differences in the behavior of 
in-plane and out-of-plane angles and their strain dependence: For all studied thicknesses of 
the STO quantum well {N = 2 — 4), the in-plane B-O-B bond angles show the strongest 
deviation from 180° in particular in the STO part, while in the LAO part the deviations 
of the bond angles are smaller and converge towards the bulk LAO value away from the 
interface. As a function of strain, the ones in the STO part {N = 2) increase with the 
lateral lattice constant, i.e. ~ 152° at ulaO) followed by ~ 155° at aNco and finally ~ 158° 
at asTO- The opposite trend is observed for the out-of-plane BOB bond angles. Here the 
strongest deviation from 180° is for osto (~ 162°), then ongo (~ 168°) and lowest for olao 
(~ 175°). Moreover, the variation in out-of-plane bond angles between the LAO and STO 
part of the SL is much weaker than the in-plane bond angles. 
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FIG. 4: In-plane (hollow squares) and out-of-plane (solid squares) Ti-0 bond lengths for the 
distinct Ti sites in (LAO)M/(STO)Ar (001) superlattices at asTo(left panel) and aLAo(right panel). 
The corresponding c/a ratio is given at the top of the graphs. 

In order to relate the structural distortions to the resulting orbital reconstruction, in Fig. 
we have also analyzed the corresponding Ti-0 bond lengths. N = 2 under tensile strain with 
c-parameter fixed to conserve the volume comprises the only case where the in-plane bond 
distance is larger than the out-of-plane Ti-0 length. This is consistent with the dxy orbital 
order at Ti^+. For all other cases, the out-of-plane Ti-0 bond distance is larger than the 
in-plane value. These differences are particularly pronounced at the Ti^+ sites, e.g. 2.10 vs. 
2.01 A at asTO and 2.13 vs. 1.97 at Ulao for = 2, respectively. A similar trend but with 











smaller absolute values is observed at the Ti^+ sites (2.02 vs. 1.96 A and 2.06 vs. 1.93 A, 
respectively). The enhanced out-of-plane to in-plane distance ratio favors the alternating 
dxz-, dyz orbital polarization. As this trend is stronger for ulaO) h is responsible for the 
stabilization of this orbital order and consequently the insulating state up to A = 3, while 
at asTO the system is already metallic for N = 3. As discussed above, the insulator-to-metal 
transition is connected with a change in orbital polarization to d^y in the interface layer and 
mixed d^z+dyz in the inner layer. Indeed the difference in out-of-plane to in-plane distance 
is reduced (ulao^ 2.06 vs. 1.93 A) and almost quenched at Ti(IF) (asTo^ 2.00 vs 1.98 A) for 
N = 4. Conversely, the inner (IF-1) layers show a larger difference (ulao: 2.10 vs. 1.92 A; 
usTO^ 2.03 vs. 1.97 A). Thus changes in the in-plane and out-of-plane bond lengths clearly 
reflect the trends in orbital polarization as a function of strain and thickness N of the STO 
quantum well. 

Chemical control: (LAO)m/(STO)7v(001) vs. 

(NGO)m/(STO)jv( 001 ) 

To explore how chemical variation of the STO counterpart influences the electronic re¬ 
construction and relaxations in the STO QW, we performed calculations on n-type 
(NGO)m/(STO)jv (001) superlattices with the in-plane lattice parameter fixed to ungO) 
which formally has the same polar discontinuity across the interface. For comparison, Figs. 

1^ and contain also the structural and electronic properties (LAO)M/(STO)Ar(001) su¬ 
perlattices at ongo- 


(a) N=2 (Onso) (b) N=3 (o^co) 



FIG. 5: Quantitative analysis of B-O-B bond angle in (OOl)-oriented (LAO)Af/(STO)iv and 
(NGO)M/(STO)Ar superlattices with the in-plane lattice parameter fixed at ongo for different 
STO thicknesses N (shaded region). {&) N — 2 STO layers, (b) A = 3 STO layers. 
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Fig. I compares the evolution of in-plane and out-of-plane B-O-B bond angles across the 
interface, for the NdGaOs/SrTiOs and LaAlOs/SrTiOs superlattices with N = 2,3. For 
N = 2 (see Fig. [^a)), the in-plane bond angle in the NdGaOs/SrTiOs heterostructure 
remains nearly constant (~ 151°) throughout the superlattice. This value is slightly smaller 
but close to the bulk value (153.2°) of orthorhombic NGO (Pbnm). For the LaAlOa/SrTiOs 
superlattice the bond angle in the STO quantum well is strongly reduced from the bulk STO 
value to ~ 155°. Moving away from the interface into the LAO part, the B-O-B bond angles 
quickly approach the LAO bulk value of 171.4°. In contrast, for both super lattices, the 
out-of-plane bond angle deviations 180° in the STO part are less pronounced and the values 
almost coincide for the NdGaOs/SrTiOs and the LaAlOs/SrTiOs superlattices in the center 
of the STO part, whereas they converge towards the NGO and LAO bulk values outside the 
STO QW. 
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FIG. 6: Side view of the relaxed structure, electron density distribution, integrated over occupied 
Ti 3d bands between F’p-l-S eV and E’p and band structure in (OOl)-oriented (NGO)m/(STO)a^ 
and (LAO)m/(STO)a 7 superlattices at ungo with N = 2. Majority and minority bands are plotted 
in blue and orange, respectively. 

Increasing the STO thickness to = 3 (see Fig. ib)), the in-plane bond angle deviations 
are signihcantly reduced in the central STO part (~ 170°), but are still quite prominent at 
the interface (~ 157°). The deviation of the out-of-plane bond angles is similar as for the 
super lattices with N = 2. 
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FIG. 7: Side view of the relaxed structure, electron density distribution, integrated over occu¬ 
pied Ti 3d bands and band structure in (OOl)-oriented (NGO)M/(STO)Ar and (LAO)M/(STO)Ar 
superlattices at ungo with N — 3. Majority and minority bands are plotted in blue and orange, 
respectively. 

The similarity in structural relaxation of the NdGaOs/SrTiOs and the LaAlOs/SrTiOs 
superlattices goes hand in hand with common electronic properties: For = 2, the ground 
state for both systems is an insulating charge-ordered phase with Ti^+ and Ti^+ sites (see 
Fig. g. Increasing the STO thickness to V = 3, the system undergoes an insulator-to- 
metal transition (see Fig. [^. The charge-ordered state is suppressed and the excess charge 
is delocalized throughout the STO slab, similar to the corresponding case at usto (see Fig. 

ic)). 

Discussion 

A central question is the origin of the unanticipated orbital reconstruction and how it relates 
to the STO QW thickness and structural distortions. Our results suggest that the electronic 
ground state and orbital polarization are intimately related to the level of electrostatic 
doping of the system. By construction the polar discontinuity at the interface introduces 
0.5e per lateral unit cell at each interface, corresponding to a carrier density of ~ 7 x 10^^ 
cm“^. N = 2 comprises an extreme case with doping of 0.5e per Ti site. With increasing 
Nsto the excess charge is redistributed throughout the STO QW, resulting in a reduction 
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of the carrier density and the effective fflling of the Ti 3d band. We consider this reduction 
of carrier density with increasing Wto as the driving force for the switching of orbital 
polarization from d^z, dyz at the Ti^"*" sites ioi N = 2 to d^y in the interface layer and mixed 
dxz, dyz in the inner layers beyond N = 3 as well as for the concomitant insulator-to-insulator 
transition. 

It is instructive to compare the LAO/STO SLs to GTO/STO systenJ^il. For the latter a 
similar critical thickness of = 3 was reported for the insulator-to-metal transition. In 
both cases the QW thickness reduction increases the effective carrier density and brings the 
system into a state where the on-site Coulomb repulsion becomes relevant, as pointed out 
by Ouellette et al^. 

Concerning the orbital reconstruction most of previous studies have reported a dominant 
dxy orbital polarizatiorP^*^. Here we find that this orbital polarization is stabilized only for 
a compressed c lattice parameter, whereas the ground state comprises a d^z, dyz ordering at 
the Ti^"*" sites. Beyond the critical thickness, a crossover to a d^y character in the interface 
layer and a mixed d^z, dyz character in the inner layers takes place. This is consistent 
with recent finding of a dominant contribution of d^z, dyz orbitals to the Fermi surface in 
LAO/STO, where doping levels beyond ~ 6.9 x 10^^ cm“^ were achieved through oxygen 
vacancy dopin^^Sl origin of orbital polarization and insulator-to-metal transition is 

further substantiated by the analysis of bond lengths (see. FigQ. The dxzi dyz orbital 
polarization at the Ti^"*" sites is associated with significant elongation of the octahedra, 
whereas the transition from insulating-to-metallic state at the critical STO thickness is 
reflected in a reduction/increase in the apical to basal bond length ratio in the IF/IF-1 
layer. 

A key structural feature are the strong tilts of the TiOe octahedra, in particular at the 
interface, in contrast to the undistorted bulk STO. In the related CTO/STO superlatticeJ^ 
this behavior was ascribed to the influence of CTO, which has bulk B-O-B angles of 155°. 
As LAO has a much larger B-O-B angle of 171.4°, we can unambiguously conclude that 
the B-O-B angle reduction in the STO QW is not induced by the STO counterpart LAO, 
but correlates with the Ti 3d band occupation and the localization of charge. This is 
furthermore confirmed by the similar structural and electronic behavior of the STO QW 
in the NdCaOs/SrTiOs (bulk NGO has B-O-B angles of 153°.) and the LaAlOs/SrTiOs 
superlattices. In particular, the strongest changes in the in-plane B-O-B angles are observed 
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at the interface at sites, whereas the changes are much smaller and even quenched in 
the inner layers or at Ti^'*' sites, as can be observed ior N = 3 and 4 in Figs, and d. 

In summary, DFT+{7 calculations, performed on (OOl)-oriented (LAO)m/(STO)jv and 
(NGO)M/(STO)iv superlattices with n-type interfaces show a rich set of orbital reconstruc¬ 
tions depending on the STO quantum well thickness and c/a ratio. For N = 2 the ground 
state is a charge ordered TF+, TF+ phase characterized by unanticipated d^zidyz orbital 
polarization and elongated Ti^^-Oe octahedra. c lattice compression leads to a Dimer Mott 
insulator (ulao) or a d^y orbital polarization (usto)) the latter case being the only one with 
compressed XF^-Oe octahedra. With increasing STO quantum well thickness an insulator- 
to-metal transition takes place. The changes in the in-plane to out-of-plane bond length 
ratios signal the switching of orbital polarization and the critical thickness for the insulator- 
to-metal transition. A central finding is the pronounced enhancement of octahedral tilts in 
the STO quantum well that are not present in the bulk and can be unambiguously associ¬ 
ated with the electrostatic doping of the polar n-type interface in these systems. Together 
with the exotic electronic states found in (lll)-oriented (LAO)m/(STO)jv SL^, the results 
demonstrate how strain and the thickness of the STO quantum well can be used to engineer 
the orbital reconstruction and insulator-to-metal transitions in oxide super lattices. 

Methods 

DFT calculations were performed on n-type (LAO)M/(STO)Ar(001) and 
(NGO)jvf/(STO)iv(001) superlattices with varying thickness M, N of the constituents, 
using the all-electron full-potential linearized augmented-plane-wave (LAPW) method, as 
implemented in the WIEN2k cod^. For the exchange-correlation functional we used the 
generalized gradient approximation (GGA)P^. Static local electronic correlations were taken 
into account in the GGA-|-[/ methocP^ with U = 5 eV, J = 0.7 eV (Ti 3d), U = 8 eV (La 
4/, Nd 4/). For oxides containing Ti^+ (3d^) values between U = 3 — 8 eV have been used 
in previous studied^®2i®. The influence of strain was explored by setting the lateral lattice 
parameter to the one of LAO (olao = 3.79A), NGO (ungo = 3.86A) or STO (ustc^ ~ 
3 . 92 A), corresponding to superlattices grown on the respective (001) substrate. Octahedral 
tilts and distortions were fully considered by relaxing atomic positions in PI symmetry 
using a \/2o x v^a lateral unit cell. Additionally the out-of-plane lattice parameter c was 
optimized for all super lattices. The effective masses of the electrons were calculated by 
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fitting the region around the Gamma point to a quadratic polynomial. 


Note added to proofs: 

Such an emergence of strong octahedral tilts associated with a transition to a Mott insulating 
phase were recently predicted also in doped bulk SrTiOj^. 
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